Complexity of spectral sequences: semiclassical approach 
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. It has been long recognized that the task of semiclassical evaluation of quantum 

spectra for the classically nonintegrable systems is fundamentally more complex than 
' for the classically integrable ones. Below it is argued that the quantum spectra of 



o 



the chaotic systems can differ among themselves by level of their complexity. 
PACS numbers: 



> ■ I. INTRODUCTION 

^ . It is well known that the structure of the spectrum of a quantum mechanical system re- 

flects the dynamical properties of its classical counterpart. The degree of classical dynamical 
regularity is reflected in the general statistical properties of the quantum spectra and in the 
nature of the analytical solution of the spectral problem. 
^ ' For the conservative systems, dynamical regularity is usually understood as the possibility 



to specify the trajectories via a set of quantities that remain constant throughout dynamical 
evolution - the integrals of motion. If each one of the d degrees of freedom of a bounded 
system corresponds to a conserved quantity, Ji, J^, then such a system is completely 
integrable and its dynamical behavior can be regarded, in suitable coordinates, as a com- 
bination of oscillations. Algebraically, this is a manifestation of the underlying symmetries 
in the system, that lead to the appearance of certain global structures in the phase space 
known as integrable tori [1], that uphold the oscillatory pattern of the trajectories. 

As it was pointed out by the EBK theory, in semiclassical regime the integrals of motion 
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assume a set of discrete values, defined via the so called quantum numbers, n^, 

h-hln, + ^), (1) 

where h is the Plank's constant and /Xj is the Maslov index [2, 3]. In view of (1), each 
integral of motion provides a uniform global map from the natural numbers into the spectral 
sequence, which leads to the analjd;ical solution to the corresponding quantum mechanical 
spectral problem. For example, if the energy is defined in terms of the motion integrals as 
E = E^^^^^^ {Ii, Id), the quantum eigenvalues of energy are 

ih, /,) ^ {n, + ^),...,h (n, + ^)) . (2) 

The integers rii can be interpreted as the number of times the wave representing the quantum 
mechanical particle winds around the basic torus cycles, so physically, the discretization of 
the spectrum in the EBK theory can be understood as a manifestation of the geometrical 
consistency between the (semiclassical) quantum waves and the dynamical trajectories. 

In contrast, if the number of motion integrals is less than the number of the degrees of 
freedom, then the dynamics is "irregular", in the sense that its trajectories do not follow 
any particular patterns in the phase space. Rather, they tend to cover as much phase space 
X as possible at a given energy, more or less uniformly, so the only relevant long term 
characteristics of the dynamics are the relative frequencies of the visitations of different 
parts of X, which can be equally well described by a certain smooth probability distribution 
Fix). 

As a result, the complexity of the semiclassical quantization task for the nonintcgrable 
systems differs substantially from the integrable ones. While the EBK quantization method 
is constructive, i.e. it is possible to define the individual levels of the quantized integrable 
systems explicitly, by specifying a particular set of quantum numbers (2), for nonintcgrable 
systems such individualized solution for the spectral problem is generally unavailable. In- 
stead, one of the main results of the semiclassical theory is the series expansion representa- 
tion for the quantum density of states, which is a global characteristics of the spectrum, the 
so-called Gutzwiller's formula: 

p{E) ^Y.^{E-En)= p{E) + Re J] B,{E)e''^^'^\ (3) 

n p 

Here Sp{E) is the action functional evaluated for the periodic orbit p, and Bp is a certain 
weight factor [2]. 
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One can immediately appreciate the difference in complexity of the series (3) compared 
to the (2). The sum (3) includes all the periodic orbits (that are isolated in a fully chaotic 
system), and so it reflects the full complexity of the classical phase space structure of a 
nonintegrable system. The essence of the Gutzwiller's formula is that the oscillating ampli- 
tudes e^^p^^^ produced by the periodic orbits, combined with the appropriate weights Bp, 
produce a constructive interference effect every time E happens to be equal to one of the 
quantum eigenvalues and cancel each other out for all other values of E, which enables one 
to transform the information contained phase space structures into the pattern of quantum 
mechanical spectral sequence [2]. 

An important aspect of (3) is that it is indeed a very general result, that can be applied 
to a great variety of systems. Moreover, certain mathematical systems or objects can be 
put into a quantum chaotic context because their constituents can be related via formula 
(3) and hence be interpreted as "quantum chaotic". One well known example of this is 
provided by the relationship between the eigenvalues of the Laplace - Beltrami operator 
on a surface of constant negative curvature and the dynamics of a point particle moving 
on it. This connection is described by the famous Selberg trace formula [4], that was 
first discovered as a number theoretic and functional analysis theorem. Apart from the 
deterministic chaotic systems, the expansion (3) can also describe the spectra of many 
classically stochastic systems, e.g. the so-called quantum graphs [5, 6] and 2D ray splitting 
biUiards [8]: as long as the stochastic dynamics generates trajectory patterns similar to 
the deterministic chaotic trajectories, they manifest themselves in the in the same physical 
context in quantum regime. 

Perhaps the most intriguing example is provided by the relationship between the zeroes 
of Riemann's zeta function and the set of the prime numbers. In [9-11] it was argued that 
the integrated spectral density, the spectral staircase of the imaginary parts of the nontrivial 
zeroes of Riemann's zeta function. En — Im(2;„), C(-^n) = 0, 

Ar(£;) = ^e (4) 

n 

can be expanded into the Gutzwiller type series (3) in which the average part of the density 
of EjiS is given by 
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FIG. 1: Spectral staircase for the imaginary parts of zeroes of Riemann's zeta function. Wiggly 
line represents the Gutzwiller's series for N{k), and the smooth line running through the staircase 
is the Weyl's average N(k). 

and the oscillating part by 

1 °° pimElnp 

TT ^ — ^ ^ — ^ mp 2 

p m=l 

where the index p runs over the set of prime numbers. Hence the prime numbers in this 
case play the role of the prime periodic orbits of lengths Lp = In p. This shows that the 
organization of the roots of the zeta function, which a is purely number theoretic object, 
is as irregular with respect to the set of prime numbers, as is the quantum spectrum of a 
classically chaotic system as expressed by the "periodic orbit expansion" [3, 12]. Despite 
these significant implications of Gutzwiller's formula, it does not produce the final solution 
of the spectral problem in the same sense as the EBK theory does for the integrable sys- 
tems. The problem is that the Gutzwiller's formula does not specify the individual energy 
or momentum eigenvalues, i.e. it does not produce a functional correspondence between 
the spectral sequence and the natural number (quantum number) sequence, although such 
correspondence actually can exist. 

In order to evaluate the individual eigenvalues En from (3), one needs to provide addi- 
tional local information, that would allow to single out the individual delta peaks in p{E). 
To have a uniform solution to the spectral problem that would be equivalent to the EBK 
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FIG. 2: Spectral staircase N{k) for the eigenvalues of momenta of quasi one dimensional quantum 
networks (quantum graphs), that are known to be stochastic in classical limit (see. Section III). 
A smooth monotone function f(k) intersects every stair-step of N{k). 

formula (2), one must be able to produce such information systematically, in a way that 
would eventually lead to an explicit functional dependence En = E{n). 

It is easy to show that in principle this not an impossible task. Indeed, consider a 
the spectral staircase (4), that monotonously increases with energy. Consider also a smooth 
monotone function f{E) that intersects every stair step of N{E) as shown on Fig. 2. It is also 
possible to chose f{E) in such a way that the average deviation of f{E) from N{E) vanishes, 
{f{E) — N{E)) = 0, so f{E) can be considered as an "average" of N{E), that is in general 
different from the Weyl's' average [16]. The intersection points En, f (j^n^ = N (j^n^ = n, 
satisfy the condition 

... ^ En-l ^ En ^ En ^ -En+1 ^ -E'n+l ^ ••• • (7) 

Since f{E) is monotone, it can be inverted, so the roots of f{En) = N (j^ri^ , can be defined 
as a function of n, En = /^^(^) = E{n). The resulting separating sequence (7) can then be 
used to find the eigenvalues En in the form of an explicit formula [20-23] , 

En = E{n) = / p{E')E'dE'. (8) 

JEin-l) 

Due to the expansion (3), one can in principle evaluate the integral (8) explicitly and obtain 
the formula for En = E[n) in an expansion form structurally similar to (3). In general, the 
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structure of the classical phase space and hence the structure of the sum (3) can change 
in very complex ways as a function of energy, which makes the evaluation of the integral 
(8) a much more complicated problem. In order to avoid these difficulties, the following 
discussion will be restricted to the so-called scaling systems, e.g. the billiards, in which case 
the expansion coefficients Bp are energy independent and the functional structure of the 
periodic orbit sum (3) is fixed. 

Such construction produces a map from the natural number sequence first into the sep- 
arating sequence — f~^{n), and then to physical spectral sequence via (8). Hence in 
fully nonintegrable system there also exists a solution to the spectral problem that produces 
the eigenvalues En as a function of n, En = E(n) based on the semiclassical expansion (3). 

Formulae (7) and (8) suggest that the possibility to solve the spectral problem depends 
on the possibility to localize the individual delta peaks within intervals In — En-, E^-i , 
loosely defined by the inequality (7). Prom such perspective, the missing part of the solution 
is the information about the structural complexity if the spectral sequence. 



II. STRUCTURAL ANALYSIS 

The task of obtaining En is in fact much more complicated than it may appear at the 
first sight. Finding a monotone function that would follow the pattern of N{k) in the detail 
required by (7) for a generic system proves to be an extremely difficult problem [16]. As a 
result, outside of a few simple systems [17-23], there are virtually no examples of explicit 
functional dependencies of the quantum energy levels, E^i on the index n, E — E{n), for 
the classically nonintegrable systems. 

The apparent impenetrability of the spectral problem for quantum chaotic systems seems 
to impose an implicit (and false) empirical dichotomy ~ either the system is integrable and 
explicitly quantizable via EBK theory, or it is nonintegrable and the spectral problem does 
not have an explicit solution. However, as it was argued above, even for a completely 
nonintegrable systems it is still possible to find the semiclassical solution to the spectral 
problem. En — E{n). Moreover, the intricacy of the quantum chaotic spectra may conceal 
a rich complexity structure, that can be studied both within the paradigm of semiclassical 
physics and outside of it. 

In itself, the mathematical problem is to describe the complexity of mapping the natural 
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numbers n — 1, 2, ... into the spectral sequence En — E{n). However, in the specific 
context of quantum chaos theory, the task is to do this within the paradigm of the semi- 
classical physics, in which the objects and the phenomena of the classical dynamics provide 
the semantics for describing quantum objects and phenomena. The question is then, how 
uniform the description of spectra may be from the point of view of algorithmic complexity? 

It is clear that a priori, the regularity of the delta peak patterns in p{k) may be different 
for different quantum chaotic systems, so the effort required for establishing a mapping from 
the natural numbers into the the sequence of E^s (and hence into the E^s) may be different, 
so the semiclassical solutions of spectral problems for the nonintegrable systems may not be 
algorithmically equivalent. 

Let us examine the regularity of the spectral sequences more closely. Let N{E) be the 
average number of the spectral points on the interval [0, E], so that N{E) — N{E) + SN{E), 
{5N{E)) — 0. For example, N{E) can be the Weyl's average, given by the volume of the 
phase space of the system for the corresponding E. This function can be used to define the 
unfolded spectral sequence 

k^:^-N{En), (9) 
distributed with uniform average density p^^\k) = 1, so that 

+ (10) 

with (^Sn^^ = 0. This unfolding operation provides a common ground for studying spectral 
sequences coming from different systems. 

Clearly, the possibility to find an analytically defined bootstrapping sequence depends on 
the magnitude of the fiuctuations 5n^ . In the simplest case, if the disorder of the original 
sequence kn^ is weak, so that it is sufficiently close to a periodic sequence, then the periodic 
points 

^^(reg) ^^ + ^(reg)^ (11) 

where ^^^^^^ is a constant, can be interlaced with it according to (7). Geometrically, this case 
corresponds to the situation when the spectral staircase N{k) of the unfolded sequence can 
be pierced by the straight line average N{k). The systems with such "almost integrable" 
spectral behavior were referred to as regular in [20-23] . The numerical values for kn^ in this 
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FIG. 3: Spectral staircase for unfolded regular spectral sequence. Every stair-step of N{k) is 
pierced by a single line. 

case can be computed explicitly given the explicit expansion of p{k), via 

= / kp^''\k)dk. (12) 

It is clear that a priori, the stair steps of a generic sequence's staircase can not all be 
pierced by a single straight line, so a single periodic sequence that bootstraps the spectrum 
does not exist. In other words, a generic sequence fci^-* is so disordered, that any sequence 
kn^ that bootstraps it cannot itself be periodic. Instead, the bootstrapping requires a 
smooth monotone curve f{k) that intersects each stair-step of the spectral staircase N{k) 
and generates a certain aperiodic sequence kn^ that is circumscribed ill IZfi ^ J cLS shown on 
Fig. 2. 

On the other hand, it is clear that if the bootstrapping function f{k) envelops the staircase 
N{k) maximally smoothly and uniformly, so that (/(/c) — N{k)) = 0, then the sequence kn ^ 
is more ordered than kn\ Since kn^ is more ordered than kn \ it may happen that kn"* 
itself can be bootstrapped by a periodic sequence, in which case kn^ will be bootstrapped 
with the periodic sequence kn^^^ = n + ^(^^s) jj^ two steps, via one auxiliary sequence kn^ . 
If however, the kn^ sequence that bootstraps kn^ must necessarily be aperiodic, then the 

(3) 

question will be whether the sequence kn can be chosen periodic, an so on. 

This immediately suggests a clear strategy of "unfolding" any sequence kii ^ using an 
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FIG. 4: Bootstrapped spectral staircase N^'^\k)s and 2 auxiliary sequences, N^^\k) and N^'^\k) 
obtained for the quantum momentum spectrum of a fully connected 4-vertex quantum network. 
The regular staircase is pierced by the line N{k). 

auxiliary set of bootstrapped seQuences kn \ kn \ kn ^\ 

^(1) < ^(2) < r,W 
Uir-l) < Uir) < uir-l) 

that starts with the original sequence kn^ and terminates when the last 
be interlaced by a periodic sequence 

n + 1 + 'j^^^s) ^ ... . At each step, the auxiliary sequences are chosen in such a way that 
the size of the fluctuations 6n^'^^ = {k[P — decreases with the increase of the index j, 
so starting from the original sequence kn \ each following sequence is closer to periodic. 
Let p^^\k) be the density of the separating points at the level j, so N^^\k) is the spectral 
staircase for the points kn\ Then the values in two bootstrapped sequences can be related 
to one another according to 

= / k'dN^J-^\k'). (14) 
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If the harmonic expansion similar to (3) is known for each sequence kn^ , 

N^){k) = N{k) + Im ^ A^^^e'^^^''^^\ (15) 

p 

with exphcitly defined A^^ and Si^\k) [5, 7, 13-15], then the integral (14) can be evaluated 
exphcitly, and so the sequence kn^ is exphcitly mapped onto kn~^\ Due to the mappings (14) 
the index n propagates through the hierarchy of bootstrapped sequences (13) and emerges 
in the 0th sequence kn \ The number of auxiliary sequences required for bootstrapping kn^ 
with kn^ therefore produces a certain complexity index for the quantum spectra that shows 
how a global index n can be consistently mapped onto the spectral sequence of any degree 
of disorder. 

There are certainly many ways in which a given sequence kn^ can be bootstrapped. The 
number of levels in the hierarchy depends on the choice of the algorithm for obtaining the 
bootstrapping sequences kn^ . Given a strategy for generating the bootstrapping sequences, 
the regularity of the sequence can be determined empirically, either by parsing through its 
individual elements or by studying the probability of occurrence of the fluctuation magni- 
tudes. In general, a separating sequence for kn^ can be written as 

4^+^) = ai^Si + (1 - (16) 

where < aH^ < 1 are arbitrary parameters. To minimize the deviation from a periodic 
sequence, one can consider the functional 

^ = E ((e^^ - - = E + (i - -nh) -ry, (17) 

n n 

where, for the fully unfolded sequences, Sn'' = ^1 + S^^^^ — Sn^^ and T = 1. Varying F with 
respect to parameters an'' under the conditions ^ ccn^ ^ 1 and kn~^^^ > k^^l\ one gets 



if ^ aii^ ^ 1, and ai^'^ = if a^'^ ^ and a^^ = 1 if a^/'^ ^ 1, so that the optimal 
bootstrapping can be explicitly defined in terms of the original sequence, which is important 
for practical studies of empirically obtained sequences, e.g. for studying the the sequences 
of numerically obtained eigenvalues of Schrodinger's operator. For crude estimates, other 
choices of etna's can be used, e.g. aH^ — 1/2 for j = 1, r. 
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for 100,000 roots of the Riemann's zeta function (found at A. Odlyzko's website, 
http : //www. dtc.umn.edu/ odlyzko/zeta .tables) and the separating sequences obtained using 
an"^ = 1/2. Since the the first sequence in which the nearest neighbor separation magnitude 
does not exceed 1 it shows that the irregularity of Riemann zeros is less than 7. If the solutions 
(18) are used, then the complexity degree of the resulting hierarchy is 2. 

The suppression of the fluctuations across the sequences k^n+i is also reflected in the sta- 
tistical properties of the deviations 6n^ . The results of numerical analysis of the spectral 
fluctuations of quantum chaotic systems show that the probability of having large fluctua- 
tions decreases with the increase of hierarchy index j (Fig 5). Although at the 0th (physical) 
level of the hierarchy the histograms of various spectral statistics show the characteristic uni- 
versal features [32, 33], the shape of the distributions at the higher levels may deviate from 
them. Numerical studies indicate that the distribution of various spectral characteristics 
at the regular level (especially for the hierarchies with high r) tends to have a Gaussian- 
like shape, which develops through the hierarchy into different physical distribution proflles. 
For example, for the nearest neighbor statistics, si^ = ^1 + 5^^}^^ — it develops into 

Wignerian-like distribution, or, for 5n^ or .^i"''' = + j /2 it develops into a Gaussian 

distribution with a larger variance, etc. A typical illustration of the appearance of the uni- 
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FIG. 6: A 6 vertex quantum graph. Dashed hue represents a periodic trajectory p with the length 



versal probability distribution profiles at the most disordered, Oth, level of the hierarchy out 
of the distributions of the the regular, rth, level is shown on Fig. 5. Clearly, the support of 
the probability distribution functions becomes progressively wider with the approach to the 
physical level of the hierarchy. At the regular level, the range of the probability distribution 
function is defined by the condition \^n^\ < 1. 



In [17-23] it was shown that this approach can be applied to the case of the quantum 
graphs - simple quasi one dimensional, scaling, classically stochastic models (Fig. 6), that 
are often used to model deterministic chaotic behavior in low dimensional dynamical systems 
[5, 6]. In [17-19] it was shown that there exist quantum graphs with different degrees of 
spectral irregularity defined in the sense of the bootstrapping hierarchy. For the simplest 
case of the regular quantum graphs [20-24], the spectrum obtained from (12) is given by 




III. COMPLEXITY OF QUANTUM GRAPH SPECTRA 




(19) 



p 
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where the coefficients Cp°^ emerge from the parameters of the Gutzwiller's series [5, 7] after 
the integration (12), Lq is the total length of the graph bonds, and L^^s are the lengths of 
the periodic orbits p on the graph. 

The expansion (19) suggests a simple physical interpretation. As mentioned above, a 
quantum graph is regular, if its spectrum is sufficiently close to the periodic sequence 

= ^- (20) 

Incidentally, (20) defines the "integrable" spectrum of a point particle in a ID box of the 
same overall size Lq. According to (19), the required corrections to (20) are given by the 
semiclassical amplitudes e'^^p \ where Sp^^ = llf^kn^ are the values of the action functionals 
evaluated for all the available periodic motions in the system with the "integrable" momen- 
tum (20). So the nth level quantum wave that is geometrically consistent with the graph's 
size, as it would be in the case infinite square well, now also circles simultaneously along all 
the possible periodic paths on the graph with the momentum (20). The resulting periodic 
orbit fiuctuation amplitudes, combined with the right weights, produce the actual momen- 
tum eigenvalues for the nonintegrable system. Hence the exact "nonintegrable" values of 
kn ^ also appear as the result of a complex interference effect similar to (3). 

This represents a generalization of the EBK quantization formula for the case of quantum 
graphs, that allows to build the individual eigenvalues kn constructively, via the explicit 
functional dependence on the index n. The effect of classical nonintegrability appears in 
(19) as a "perturbation" to the integrable background pattern, if /c„ — ^ < In this 
aspect, the situation is reminiscent of the effect of persistence of the integrability structures 
in the phase space under small nonintegrable perturbations described by the KAM theory 

[!]• 

For the irregular graphs with several level spectral hierarchy, the relationship (19) is 
repeated for each transition between the hierarchy levels, 

k^) = 4^+1) + r7(^+^) + Im J] C^)e^^'^^'' , (21) 

p 

where r}^^^^") is a bounded function of the fiuctuations on the level, ■q'^^^^'> — 7^0+^) (5(-^+^)) 
[25-27]. Also in general, the harmonic expansion (21) is different from the periodic orbit 
expansion of (19), C^^^ ^ C^^ and L^J^ ^ L^J^ [25-27]. Hence the jth level oscillations 
transmit the discrete momentum values from the {j -|- l)th to the jth level of the hierarchy. 
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As a result, the "integrable spectrum" kn' that exphcitly carries the quantum number n is 
adjusted r times according to the set of equations (21) until it is transmitted from the regular 
to the 0th level of the hierarchy. In [17-19] it was sown that every graph is characterized by 
a finite irregularity degree. 

It is important that a quantum graph of given topology can have different degrees of 
spectral irregularity depending on the bond length and other graph parameters [17-19]. 
Since the geometric complexity of the periodic orbit set is defined by the topology of the 
graph, it means that spectral irregularity, as a complexity measure, is not a trivial reflection 
of the phase space complexity of the underlying classical system, and provides a separate 
characterization of the complexity of quantum spectra. 

The expansions (19) and (21) can be used to describe analytically the development of the 
universal probability distributions [25-27], illustrated on Fig. 5. 

IV. DISCUSSION 

The task of quantifying the complexity of the map between the natural numbers and a 
given sequence x„ = x{n) is very general. This problem was recently considered in [30, 31] 
for the case of finite sequences over finite alphabets, which revealed a remarkably complex 
organizational structure of these mathematical objects. The complexity organization scheme 
developed in [30, 31] is generated by the hnear difference operator 

Ax{n) = x{n) - x{n + 1) = x'{n), (22) 

that maps one sequence into another. The motivation for using the operator A is that in 
certain simple cases it restricts the complexity of the symbolic sequences, i.e. it produces 
more ordered sequences out of less ordered ones. For example, the constant sequences 
Xn = Xn^ — const, are mapped into = Axn^ — 0, which is the "simplest" constant 
sequence. If x„ is linear, x„ = Xn^ — an + b, then its image x' — Axn^ — Xn^ is a constant 
sequence, which is simpler than linear, so that A^x^^^ = 0, and so on. In general, the 
functions for which A'x = for some / are the polynomials of degree m < I, x = p^"*) e Pi. 
Intuitively, the higher is the degree of the polynomial, the more complex is the mapping 
Xn = x{n). 

For the case of finite sequences, one necessarily runs not only into the polynomials. 
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but also into more complex "exponential" functions. By definition, the function for which 
A'x(n) = x{ri) is an exponential polynomial e^^^ e Ei of the order q that divides I. In 
[30, 31] it was shown that any function / can be characterized by a polynomial degree 
dcg(/) = deg (p*^™)) = m and an exponential order ord(/) = ord (e'-^-') = q, which allows to 
formalize the organization of complexity of the sequences. 

By definition [30, 31], a function Xn is more complex than if ord(a;„) > ord(x^). If 
two sequences have the same order, ord {x'^) — ord {x'^) then x„ is more complex than x'^ if 
deg{xn) > deg(x^). This definition of the complexity scales is similar to the organization 
of the growth rates of the exponential polynomials as described e.g. in [34] . Empirical (e.g. 
numerical) studies of simple examples, e.g. of the binary sequences, show that it captures the 
intuitive idea that, e.g., 010101010101 is simpler than 010001101111, and helps to establish 
a number of beautiful relationships and a surprisingly rich complexity structure. 

It is clear however, that the case of finite discrete sequences over finite alphabets is 
simpler than the case of the infinite sequences with real valued elements, as in the case of 
the spectral sequences. Although in some cases it is possible to study the complexity of 
discretized sequences [29], for spectral sequences it is not clear a priori which discretization 
scheme should be used. 

Using the bootstrapping hierarchy approach leads to a natural scale of complexity for 
the spectral sequences. The bootstrapping transformations suggested by the use of the 
Gutzwiller's trace formula, also produce more ordered sequences out of less ordered ones, 
and generate a finite complexity degree for a number of systems, such as spectral sequences 
of quantum graphs or the nontrivial zeroes of the Riemann's zeta function. This degree is 
analogous to the polynomial degree of the finite sequences as defined in [30, 31]. 

However, an important feature of the finite sequence complexity structures established 
in [30, 31] via (22) is that they do not generalize in any trivial way when the length of a 
sequence is increased. On the one hand, in many applications, a given A'"-letter long sequence 
over an M-letter alphabet may appear as an approximation to a longer (e.g. infinite) A^i- 
letter long sequence, A^i > iV, defined over a larger Mi-letter alphabet. Mi > M, and the 
task is to characterize the complexity of the entire sequence. Hence it is natural to use 
complexity measures that are stable with respect to such completions of shorter sequences 
by the longer ones. In the context of studying spectral sequences, complexity index defined 
for a sufficiently long list of elements should stabihze with the increase of sequence's length. 



16 



so that the spectrum as a whole is characterized by a single coherent complexity degree, 
that can be deduced from sufficiently long finite approximations. 

In view of this, it is particularly significant that the degree produced by the bootstrapping 
method proposed above allows a stable characterization of the complexity of the whole 
spectral sequence. Numerical analysis of Riemann's zeroes and of the spectra of quantum 
graphs of different topologies shows that the regularity degree obtained on the basis of a few 
hundred levels remains the same when a much larger (10^ — 10^) set of levels is considered. 

It is also physically relevant that the expansion (3) is typically derived with a semiclassical 
accuracy, so the locations of the peaks generated by the sum (3) corresponds to the actual 
En values only approximately. Moreover, even in cases when Gutzwillcr's formula is exact, 
as in the case of the quantum graphs [5, 7] and a few other systems [2, 4, 10, 35] the sum (3) 
cannot in general be computed exactly, because only a finite number of the periodic orbits 
may be known and their characteristics (e.g. the factor coefficients Bp and the actions Sp) 
are described with a finite accuracy. Hence the harmonic expansions on each level of the 
hierarchy provide only a "fuzzy" description of the spectrum. It is therefore important that 
the regularity degree obtained via the bootstrapping method is also stable with respect to 
the using finite approximations to Gutzwillcr's sum, so the semiclassical description gives a 
correct estimate of the exact hierarchy index of complexity of quantum spectra. 

The work was supported in part by the Sloan and Swartz Foundations. 
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